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Einstein’s general theory of relativity from 1915' remains the most successful 
description of gravitation. From the 1919 solar eclipse’ to the observation of 
gravitational waves’, the theory has passed many crucial experimental tests. However, 
the evolving concepts of dark matter and dark energy illustrate that there is much to 
be learned about the gravitating content of the universe. Singularities in the general 
theory of relativity and the lack of a quantum theory of gravity suggest that our 
picture is incomplete. It is thus prudent to explore gravity in exotic physical systems. 
Antimatter was unknown to Einstein in 1915. Dirac’s theory* appeared in 1928; the 
positron was observed? in 1932. There has since been much speculation about gravity 
and antimatter. The theoretical consensus is that any laboratory mass must be 
attracted‘ by the Earth, although some authors have considered the cosmological 
consequences if antimatter should be repelled by matter’ ”°. In the general theory 

of relativity, the weak equivalence principle (WEP) requires that all masses react 
identically to gravity, independent of their internal structure. Here we show that 
antihydrogen atoms, released from magnetic confinement in the ALPHA-g apparatus, 
behave in a way consistent with gravitational attraction to the Earth. Repulsive 
‘antigravity’ is ruled out in this case. This experiment paves the way for precision 
studies of the magnitude of the gravitational acceleration between anti-atoms and 
the Earth to test the WEP. 


The weak equivalence principle (WEP) has recently been tested for 
matter in Earth’s orbit” with a precision of order 10. Antimatter has 
hitherto resisted direct ballistic tests of the WEP due to the lack of a 
stable, electrically neutral, test particle. Electromagnetic forces on 
charged antiparticles make direct measurements in the Earth’s gravi- 
tational field extremely challenging”. The gravitational force ona pro- 
ton at the Earth’s surface is equivalent to that from an electric field of 


about 10” V m”. The situation with magnetic fields is even more dire: 
acryogenic antiproton” at 10 K would experience gravity-level forces 
ina magnetic field of order 10° T. Controlling stray fields to this level 
to unmask gravity is daunting. Experiments have, however, shown that 
confined, oscillating, charged antimatter particles behave as expected 
when considered as clocks“ in a gravitational field. The abilities to 
produce” and confine" antihydrogen now allow us to employ stable, 


‘Department of Physics and Astronomy, Aarhus University, Aarhus, Denmark. Department of Physics, Faculty of Science and Engineering, Swansea University, Swansea, UK. School of Physics 


and Astronomy, University of Manchester, Manchester, UK. “Cockcroft Institute, Sci-Tech Daresbury, Warrington, UK. "University of Brescia, Brescia and INFN Pavia, Pavia, Italy. TRIUMF, Vancouver, 


British Columbia, Canada. ‘Instituto de Fisica, Universidade Federal do Rio de Janeiro, Rio de Janeiro, Brazil. Department of Physics, University of California at Berkeley, Berkeley, CA, USA. 


°Department of Physics and Astronomy, University of British Columbia, Vancouver, British Columbia, Canada. Accelerator and Technology Sector, CERN, Geneva, Switzerland. "Department 


of Physics and Astronomy, York University, Toronto, Ontario, Canada. Department of Physics and Astronomy, University of Calgary, Calgary, Alberta, Canada. “Department of Physics, Simon 


Fraser University, Burnaby, British Columbia, Canada. “Van Swinderen Institute for Particle Physics and Gravity, University of Groningen, Groningen, The Netherlands. "Department of Physics, 


Stockholm University, Stockholm, Sweden. ‘Department of Physics, British Columbia Institute of Technology, Burnaby, British Columbia, Canada. "Department of Chemistry, University of 


British Columbia, Vancouver, British Columbia, Canada. “Institute for Particle Physics and Astrophysics, ETH, Zurich, Switzerland. Experimental Physics Department, CERN, Geneva, Switzerland. 


2°Department of Physics and Astronomy, Purdue University, West Lafayette, IN, USA. Accelerator Systems Department, CERN, Geneva, Switzerland. ?Soreq NRC, Yavne, Israel. Department of 


Physics, Ben Gurion University, Beer Sheva, Israel. “INFN Pisa, Pisa, Italy. School of Mathematical and Physical Sciences, University of Sussex, Brighton, UK. Physics Department, Marquette 
University, Milwaukee, WI, USA. “e-mail: william.bertsche@cern.ch; joel@physics.berkeley.edu; jeffrey.hangst@cern.ch 


716 | Nature | Vol 621 | 28 September 2023 


Analysis coil 


Mini-octupole 


A LOc 


Transfer coil 


Upper MCP i 
fa ay g 
eC af i Mirror coil G 
Ssa mmm Al mn eaae A (with bias) 
Main solenoid OcB 
Barrel . $ 
scintillator aus 
Mirror coil E 
Time projection 
chamber Mirror coil D 
Mirror coil C 
Penning ti vi 
enning trap ‘ i 
PEA Mirror coil B 
Liquid helium g , 
Soluma Mirror coil A 
Background 
aer || Be Y coil 
Lower MCP and ir Capture 
electron source solenoid 
Penning 
h trap 
H o electrodes 


— ]- 
Antiprotons Positrons 


Fig. 1 | ALPHA-g apparatus. a, Cross section of the ALPHA-g apparatus. The 

full device comprises three antihydrogen trapping regions; only the bottom 
one is employed here. The MCP detectors are used to image charged particles 
(e, e*,p) extracted from the Penning traps for diagnostic purposes. b, Expanded 
view of the bottom antihydrogen trap (the dashed rectangle in a) illustrating 
the Penning trap for antihydrogen production and the superconducting 

coils that form the neutral atom trap. The on-axis, axial field profile at full 
current is shown on the right. Note that the rTPC, the barrel scintillator andthe 


neutral anti-atoms in dynamic experiments where gravity should play 
a role. Early considerations””° and a more recent proof-of-principle 
experiment” in 2013 illustrated this potential. We describe here the ini- 
tial results of a purpose-built experiment designed to observe the direc- 
tion and the magnitude of the gravitational force on neutral antimatter. 


Antihydrogen and ALPHA-g 


Trapping and accumulation” of antihydrogen are now routine, with up 
to several thousand atoms having been simultaneously stored in the 
ALPHA-2 device”. To date, all of the measurements of the properties 
of antihydrogen™ ” have been performed in ALPHA magnetic traps. In 
2018, the ALPHA-g machine—a vertically oriented antihydrogen trap 
designed to study gravitation—was constructed. The experimental 
strategy is conceptually simple: trap and accumulate atoms of anti- 
hydrogen; slowly release them by opening the top and bottom barrier 
potentials of the vertical trap; and try to discern any influence of gravity 
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main solenoid are not drawn to scale here; see Fig. 1a for a scaled image. 

The mirror coils B-F, the analysis coil, the mini-octupole, the transfer coil and 
the background coil are not utilized here. The capture solenoid is used for 
charged particle transfer and manipulations and is de-energized for gravity 
measurements. The LOc coils (dark blue in the figure) extend past the trapping 
region used here and constitute part of two additional antihydrogen traps 
intended for future use. 


on their motion when they escape and annihilate on the material walls of 
the apparatus. The trapped anti-atoms are not created at rest but havea 
distribution of kinetic energies consistent with the trap depth of about 
0.5 K (we employ temperature-equivalent energy units). Gravity is 
expected to be manifested as a difference in the number of annihilation 
events from anti-atoms escaping via the top or the bottom of the trap. 

The experimental layout is shown in Fig. 1. Antiprotons from the 
CERN Antiproton Decelerator”’ and the ELENA (Extra Low ENergy 
Antiproton)” ring are first caught in a separate, high voltage Penning 
trap in a 3 T solenoid magnet (not shown). ELENA typically delivers 
7.5 x 10° antiprotons at 100 keV every 120 s. About 5 x 10° of these are 
dynamically captured. After being cooled by co-trapped electrons, 
antiprotons are injected into ALPHA-g and dynamically re-trapped. 
A superconducting solenoid provides the background field of 1 T for 
confining the charged particles. Positrons froma Surko-type accumula- 
tor” are also injected into ALPHA-g and re-trapped; there are typically 
3 x 10° available for each mixing cycle with antiprotons. The beamline 
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Fig. 2 | Illustrations of the magnetic bias. a, Expanded view of the end-of-ramp 
mirror coil peak regions for a bias of -1g (note the discontinuous abscissa). The 
square points represent offline ECR measurements carried out to determine 

the field profile and to find the peak field location. The points with red circles 
indicate the axial locations at which ECR measurements were made at the 


for guiding the bunches of positrons and antiprotons into ALPHA-gis 
described elsewhere”. Following manipulations to control their size 
and density**, the positron plasmas are mixed with antiproton plasmas 
in a region (electrodes B23 to B35 in Fig. 1) situated within the super- 
conducting antihydrogen trap. The anti-atom trap comprises octupole 
magnets for transverse confinement and two solenoidal ‘mirror coils’ 
(Aand Gin Fig. 1) for axial (vertical) confinement. Antihydrogen atoms 
produced with sufficiently low kinetic energy can be trapped due to 
the -p-B interaction of their magnetic moments with the external fields. 
For the field strengths in ALPHA-g, the anti-atoms are spin-polarized, 
and the scalar magnitude of the magnetic field determines the trapping 
potential. The entire production and trapping region is cooled to near 
4 K bythe liquid helium bath for the trap magnets. ALPHA-g currently 
traps a few antihydrogen atoms per mixing cycle, but antihydrogen 
atoms can be accumulated” over many cycles from ELENA. We refer 
to this process as ‘stacking’. The atom trapping volume is nominally a 
vertical cylinder of 4.4 cm diameter and 25.6 cm height. 


The effect of gravity 


The experimental protocol was to stack antihydrogen atoms, then 
release them by ramping down the current in the two mirror coils simul- 
taneously over 20 s. The anti-atoms could escape either to the top of 
the trap (through mirror G) or the bottom (through mirror A) and sub- 
sequently annihilate on the walls of the apparatus (Fig. 1). The annihila- 
tions and their positions (vertices) could be detected and reconstructed 
using the ALPHA-g radial time projection chamber (rTPC) detector 


718 | Nature | Vol 621 | 28 September 2023 


beginning and end of the mirror coil ramp-down for each gravity trial. 

b, Calculated on-axis final well shapes (after ramp-down) for the positive 
bias trials. The features at |z| > 20 cm are due to the OcB (Fig. 1) end turn 
windings. The vertical dashed lines represent the physical axial midpoints 
of mirrors AandG. 


(Fig. land Methods). A coaxial, barrel-shaped scintillator detector was 
also used for event selection (Fig. 1and Methods). 

Numerical simulations of atom trajectories (Methods) indicate that 
if hydrogen atoms were trapped and gradually released from a verti- 
cally symmetric trap (that is, the on-axis magnetic field maxima are 
equal; B, = Bc) under ALPHA-g conditions, about 80% of them would 
exit through the bottom, the asymmetry being due to the downward 
force of gravity. The goal of the current experiment was to test this 
behaviour for antihydrogen. Vertical gradients in the magnetic field 
magnitude can obviously mimic the effect of gravity. Quantitatively, 
the local acceleration of gravity g, which is about 9.81 ms”, is equiva- 
lent toa vertical magnetic field gradient of 1.77 x 10°T m”acting ona 
hydrogen atom in the ground state. The peaks in the mirror coil axial 
field strength are separated by 25.6 cmat full current, soa field differ- 
ence of 4.53 x 10 T between these points would mimic gravity. This 
consideration sets the scale for the required degree of magnetic field 
control for this experiment, but it also allows us to refine the simple, 
symmetric release procedure to more systematically probe gravity. In 
particular, it is possible to either counteract or supplement gravity by 
introducing a differential current to one of the mirror coils. 

We first consider a simplified, one-dimensional on-axis model. As 
the mirror fields are ramped down, a particular anti-atom will escape 
when its axial kinetic energy exceeds the combined gravitational 
and magnetic potential at the peak axial field position of one of the 
mirror coils. Thus, one could balance the effect of gravity on matter 
by imposing a field difference (B; - B,) of about -4.53 x 10 *T between 
the mirror field peaks (Fig. 2a). Maintaining this difference during the 
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Fig. 3 | Escape histograms. The raw event z-distributions are displayed as 
histograms for each of the bias values, including the +10g calibration runs. 
These are uncorrected for background or detector relative efficiency. The time 
window represented here is 10 s to 20 s of the magnet ramp-down. The z-cut 
regions are indicated by the solid, diagonal lines. Explicitly, the acceptance 
regions in zare [-32.8, -12.8] and [12.8, 32.8] cm for the ‘down’ and ‘up’ regions, 
respectively. 


ramp-down would in principle result in half of the atoms escaping in 
each direction. Note that this incremental field is very small compared 
to the size of the initial peak end field, whichis about 1.74 T. The mirror 
coils A and G were connected in series, and a bipolar current supply 
connected only to mirror G could provide a field increment or decre- 
ment (Methods). We emphasize that a magnetic gradient is not applied 
uniformly over the length of the trap. The local field geometry in the 
region of each mirror coil determines which particles can escape axially. 


The release experiment 


Inanticipation of future precision experiments, the octupole fields in 
ALPHA-g can be generated by three distinct coils. Two of these, which 
we designate long octupole (LOc) and bottom octupole (OcB), are 
employed here (Fig. 1). The OcB magnet (made up of six wound current 
layers) spans the axial trapping region employed in the current experi- 
ment. The LOc magnet comprises two layers of windings and extends 
over 1.5 m of the apparatus, covering two additional antihydrogen 
trapping regions not utilized here. For trapping and stacking, both 
octupole magnets are energized to about 830 A. At the completion 
of stacking, the LOc magnet is ramped down in 1s, thereby eliminat- 
ing the transverse confinement field above mirror G (Fig. 1). This step 
releases some of the more transversely energetic atoms - about half of 
the stacked sample. By counting the resulting annihilations, we obtain 
an indication of the total number of atoms that have been stacked. 

The actual experiment involved many trials of antihydrogen accu- 
mulation and release for various magnetic ‘bias’ levels. We define the 
imposed bias as: 


H,(Bg ~ By) 
My(Z¢- Za) 


where H,is the Bohr magneton, (B, - B,) is the difference between the 
on-axis field maxima under the two mirror coils, my is the hydrogen 
gravitational mass and (zç - Z,) is the height difference between the 
positions of the on-axis field maxima. It is convenient to express the 
bias relative to g. Thus, inthe one-dimensional model, a magnetic bias 
of -1gwould effectively balance the downwards gravitational force for 
hydrogen. Having assumed no a priori direction or magnitude for the 


gravitational force on antihydrogen, we investigated nominal bias 
values of +3g, +2g, +1.5g, +1g, +0.5g and Og. Figure 2b illustrates the 
positive bias fields (B, > B,), which would encourage antihydrogen 
atoms to exit at the bottom. 

We typically accumulated anti-atoms for 50 stacks in roughly four 
hours, resulting in about 100 atoms trapped. For each trial, after the 
conclusion of stacking and the LOc ramp-down, the on-axis field 
magnitude at one axial location under each mirror coil (Fig. 2a) was 
measured using the technique of electron cyclotron resonance (ECR) 
(Methods). The ECR measurement was made at approximately 130 s 
after the LOc ramp-down. The mirror coil current ramp-downs hap- 
pened next and were linear over 20 s. The smaller of the two mirror 
fields was not ramped all the way down to the level of the bottom of 
the confinement well but stopped at about 5 x 10° T above this level. 
This was to ensure that the released atoms possessed enough energy 
to overcome the small axial field bumps that arise from the end wind- 
ings of the OcB magnet (Fig. 2b). At approximately 96 s after the mirror 
ramp-down, the ECR measurements were repeated to characterize the 
final axial well (Methods). 

Various bias values were interleaved during the data-taking 
period, which lasted about 30 days. We emphasize that the integer or 
half-integer bias values identified above are just labels for the trials and 
refer to the programmed on-axis field maxima; neither is the bias per- 
fectly constant during the ramp-down, nor does the one-dimensional 
model completely characterize the three-dimensional experiment. 
Trials for a given bias were repeated six or seven times, depending on 
the total number of events detected. The raw results (no background 
subtraction or detector efficiency correction) are presented as axial 
annihilation distributions in Fig. 3. For further analysis, we exclude 
events whose z-position lies between the physical mirror centres, or 
more than 0.2 m outside the physical mirror centres, as indicated in 
Fig. 3. This ‘z-cut’ was chosen by conducting a separate set of experi- 
ments in which we attempted to release trapped antihydrogen atoms 
to only the top or the bottom of the trap by applying a bias of -10g or 
+10g, respectively. The +10g trials also help to determine the relative 
efficiency of the rTPC detector for the up and down escape regions 
(Methods). The efficiency determination uses the number of atoms 
detected in the LOc ramp-down as a normalization. The plotted event 
distributions were also subject to a ‘time cut’: events are accepted from 
10 to 20 s of the ramp-down, as we found that the number of atoms 
emerging before 10 s is negligible (Fig. 4). 
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Fig. 4 | Time structure of the annihilation events from escaped antihydrogen. 
The number of detected events (left ordinate) is plotted as a function of time as 
the magnets are ramped down. This figure represents the sum of the seven trials 
having bias Og. The dashed line (right ordinate) illustrates the calculated axial 

well depth during the magnet ramp-down. The excluded events fail the time cut. 
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Table 1| Results of the release trials 


Nominal bias (g) Number of N., (events) Na (events) Events during 


trials LOc ramp-down 

-3.0 7 151.7 16.5 199.2 
-2.0 7 128.7 33.5 168.2 
-1.5 6 128.9 577 192.0 
-1.0 7 69.7 62.5 183.2 
-0.5 7 55.7 67.5 201.2 
(0) 7 36.7 94.5 144.2 
0.5 7 36.7 124.5 177.2 

1.0 7 17.7 119.5 185.2 
1.5 6 13.9 180.7 234.0 
2.0 7 67 163.5 228.2 
3.0 7 77 147.5 199.2 
-10.0 6 142.9 0.7 169.0 
10.0 6 -0.1 185.7 213.0 


The number of events for anti-atoms escaping either up or down is tabulated for each bias 
series. These events occur in the time window 10-20s during the ramp-down and lie within 
the z-regions illustrated in Fig. 3. Also shown is the number of events due to antihydrogen 
atoms that escape when the long octupole magnet is ramped down. All values are corrected 
for the expected cosmic ray background. Counting uncertainties are not listed but are used 
in the global determination of Pan in Fig. 5. The background per trial was 0.18+0.01 events in 
the top region and 0.21+0.01 events in the bottom region. The background per trial for the 
LOc ramp-down window (duration 13.1s) was 0.83+0.02 events. The +10g entries are for the 
calibration trials (see text). 


The essential cumulative result for each bias can be represented by 
two numbers, N,, and Nar: the number of particles escaping upwards 
or downwards. These are listed in Table 1. The techniques used to maxi- 
mize the signal and suppress the background are described in Methods. 
The background rates are listed in the Table 1 notes. 


The escape curve 


In Fig. 5 we plot the probability for an antihydrogen atom to escape 
downwards (P,,,) as a function of the applied bias. The probabilities 
and their credible intervals were obtained from the raw event counts 
by using standard statistical techniques (Methods). The biases plot- 
ted here are derived values, as the magnetic field difference (on axis) 
between the upper and lower barriers remains only approximately 


constant as the current is decreased. This is due to small asymmetries 
in the background field, the construction of the mirror coils and the 
ramp-induced persistent currents in the superconductors (Methods). 
We also observe that these currents decay after the end of the ramp 
(Extended Data Fig. 6), affecting the final-well ECR measurements. To 
account for these effects, we use a measurement-based magnetic field 
model (Methods) to calculate the bias during the ramp. We can then 
assign to each annihilation event the calculated bias for the time at 
which that particular anti-atom escaped the trap (Fig. 4 and Extended 
Data Fig. 8). Finally, we average the biases for all of the events that pass 
our selection criteria (or ‘cuts’) to arrive at the plotted bias value for 
the collection of trials sharing the same magnetic field configuration. 
The uncertainties in the bias determination are of order 0.1g and are 
described in detail in Methods. 

Qualitatively, the experimental data in Fig. 5 exhibit the behaviour 
characteristic of gravitational attraction between antihydrogen and 
the Earth. At a bias of about +3g(-3g) the anti-atoms exit predomi- 
nantly at the bottom(top) of the trap, as the magnetic imbalance is 
significantly larger than lg. The fraction exiting through the bottom 
increases monotonically as the bias increases from -3g to +3g. The 
balance point (Pin = 0.5) is close to -1g, as naively expected from the 
simplified one-dimensional argument presented above. 

To gain more quantitative insight into the results (and origi- 
nally to inform the design of the experiment) we rely on extensive 
numerical simulations (Methods) of the trajectories of antihydrogen 
atoms trapped and then released. The numerical model features a 
three-dimensional magnetic field map based on both the as-built 
superconducting magnet wire model and the measured fields from 
ECR or a magnetron frequency measurement technique (Methods). 
The actual currents measured during the experimental sequence 
are used for the simulation. This is the same magnetic field model 
used to derive the plotted biases above, so the simulation describes a 
three-dimensional system that is consistent with our best experimental 
measurements—both static and dynamic—of on-axis field strengths. 
The ECR measurements taken during the trials have been supple- 
mented by extensive offline studies using both ECR and the magnetron 
method (Methods). The simulated release results are plotted with 
the data in Fig. 5, both for attractive (normal) gravity and, by way of 
comparison, for ‘no’ gravity and for ‘repulsive’ gravity. 

The agreement between the shape of the measured data and that of 
the simulation is visually compelling. To extract a value for the local 
acceleration from our dataset, we have compared the data to a set of 
simulations that presume values for antihydrogen’s gravitational 
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Fig. 5 | Escape curve and simulations. The derived P,,, values are plotted versus 
bias for the experimental data and for simulations of the experiment for three 
values of the gravitational accelerationa,:1g (normal gravity, orange), Og (no 
gravity, green) and -1g (repulsive gravity, violet). See the text for the definitions 
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of the uncertainties. The right ordinate is the down-up asymmetry A = 2P,,— 1. 
The confidence intervals on the no- and repulsive gravity simulations are 
comparable to those for the normal gravity simulation and have been omitted 
for clarity. 


Table 2 | Uncertainties in the bias determination 


Table 3 | Uncertainties in the determination of a, 


Uncertainty Magnitude (g) Uncertainty Magnitude (g) 
ECR spectrum width 0.07 Statisticaland Finite data size 0.06 
Repeatability of (Bg - By) 0.014 systematic Calibration of the detector efficiencies inthe 0.12 

Peak field size and z-location fit 0.009 upand downregions 

Field decay asymmetry (A to G) after ramp 0.02 - - Other aus Soules: — - gol 

Bias variation intimne 0.02 SH pines the magnetic fields (on-axis 0.16 

Field modelling Be Antihydrogen initial energy distribution 0.03 


Summary of the uncertainties in the derived bias values, expressed in units of the local 
acceleration of gravity for matter (9.81ms~”). See Methods for definitions and details. 


acceleration that differ from 1g (Extended Data Fig. 1). Generally speak- 
ing, the simulated curves have the same shape and are shifted along 
the bias axis. From a likelihood analysis (Methods) onthe experimen- 
tal data, we find that the local gravitational acceleration of antihydro- 
gen is directed towards the Earth and has magnitude ag = (0.75 + 0.13 
(statistical + systematic) + 0.16 (simulation))g, where g = 9.81 m s”. 
Within the stated errors, this value is consistent with a downward 
gravitational acceleration of 1g for antihydrogen. 


Classification of uncertainties 


Broadly speaking, we characterize three different types of uncertainty. 
The uncertainties regarding magnetic field measurement and model- 
ling affect the derived bias values and are listed in Table 2 and described 
in Methods. These are reflected in the horizontal error bars on the bias 
values in Fig. 5. Statistical and systematic uncertainties regarding event 
detection, suchas counting statistics, backgrounds and detector effi- 
ciencies, are listed in Table 3. These are manifested as vertical error 
bars in the P,, values in Fig. 5. Finally, an estimated uncertainty band 
(orange band in Fig. 5) is associated with the simulation. This includes 
the potential impact of various unmeasured quantities, such as magnet 
winding misalignments, off-axis persistent magnetic fields, and uncer- 
tainty in the energy distributions (longitudinal and transverse) of the 
trapped antihydrogen atoms. All of the above are used to extract the 
uncertainties in the quoted value of aṣ. Our goal here is not to make a 
precision determination of the magnitude of aș, but to identify the 
statistical sensitivities and systematic effects that will be important 
for future measurements. 

As a cross check, we conducted trials in which we used a 130 s ramp- 
down time, for biases of Og, -1g and —2g. Within the calculated uncer- 
tainties, the results were consistent with the 20 s data and with the 
appropriate simulation (Extended Data Fig. 2). 

We also observe that some atoms are released after the end of the 
20 s ramp (Fig. 4 and Extended Data Fig. 3). This is potentially due to 
long-time-scale mixing’ between the transverse and longitudinal 
motions of the atoms, but this has not yet been investigated in detail. 
The gravitational behaviour of these atoms appears to be consis- 
tent with the 20 s ramp-down sample (Extended Data Fig. 3), but the 
detailed systematic measurements to confirm this have not yet been 
performed. 


Possible complicating effects 


We have considered other effects that could mimic a gravitational force 
or add significant uncertainty, and we can rule them out due to their 
negligible magnitudes. We have earlier determined an experimental 
limit for the antihydrogen charge” to be less than about 10°C. Thus, 
a1V potential change would have the same effect as a10°T change 
in magnetic field. The trap electrodes are maintained at their com- 
mon ground to within +10 mV after stacking is completed, so even 
the extremely unlikely presence of the maximal non-zero charge on 


Summary of the uncertainties involved in the determination of the gravitational acceleration 
ag The uncertainties are one standard deviation and are expressed in units of the local 
acceleration of gravity for matter (9.81ms~). See Methods for the details. 


antihydrogen would play no role here. Concerning the size of the 
magnetic dipole moment of antihydrogen, we earlier measured the 
microwave transition” within the hyperfine-split ground state at 
approximately 1 T with an absolute accuracy corresponding to 0.3 mT. 
Since the positron magnetic dipole moment mainly determines the 
transition frequency, this corresponds to an uncertainty of the mag- 
netic dipole moment of less than 1 part per thousand in antihydrogen, 
leading to anegligible contribution to the error budget here. 

The measured masses and charges of the positron and antiproton’? 
can, in the absence of new physics, be used to constrain the polariz- 
ability of an antihydrogen atom in the ground state to approximately 
that of the hydrogen ground state”: 7.4 x 10C? (J m)”. Thus, a change 
in electric field of 100 V m™ would have an effect equivalent toa change 
in magnetic field of less than 10“ T. Finally, antihydrogen atoms may 
change their velocity due to collisions with background gas during 
the ramp-down. From the measured antiproton storage lifetime of 
4,000 sin the trap, we estimate the density of background gas to be 
approximately 2 x 10°cm”. Using this value together with the calcu- 
lated cross sections*®, the probability for a collision during the 20 s 
(130 s) ramp-down is less than 0.5% (3%). 


Conclusion 


We have searched for evidence of the effect of gravity on the motion of 
particles of neutral antimatter. The best fit to our measurements yields 
a value of (0.75 + 0.13 (statistical + systematic) + 0.16 (simulation)) g 
for the local acceleration of antimatter towards the Earth. We conclude 
that the dynamic behaviour of antihydrogen atoms is consistent with 
the existence of an attractive gravitational force between these atoms 
and the Earth. From the asymptotic form of the distribution of the 
likelihood ratio as a function of the presumed acceleration, we esti- 
mate a probability of 2.9 x 10“ that a result, at least as extreme as that 
observed here, could occur under the assumption that gravity does 
not act on antihydrogen. The probability that our data are consistent 
with the repulsive gravity simulation is so small as to be quantitatively 
meaningless (less than 10). Consequently, we can rule out the exis- 
tence of repulsive gravity of magnitude 1g between the Earth and anti- 
matter. The results are thus far in conformity with the predictions of 
General Relativity. Our results do not support cosmological models 
relying on repulsive matter—antimatter gravitation. 


Future perspectives 


This experiment marks the beginning of detailed, direct inquiries into 
the gravitational nature of antimatter. Having determined the sign and 
approximate magnitude of the acceleration, our next challenge is to 
extend the method to measure the magnitude as precisely as possible, 
to provide a more stringent test of the WEP. Colder atoms will obviously 
allow for more sensitive measurements, and our simulations indicate 
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that using colder antihydrogen atoms will in general steepen the transi- 
tion region of the escape curve and allow for higher precision. Our recent 
demonstration of laser cooling of trapped antihydrogen” is a promising 
developmentin this direction. Additionally, our future measurements 
willincorporate adiabatic expansion cooling of trapped antihydrogen*. 
In addition to future measurements in ALPHA-g, alternative approaches 
are being pursued by the GBAR” and AEgIS* collaborations at CERN. 

The dependence on simulations is not a concern at the current level 
of precision, but supplementary experiments to benchmark and refine 
the simulations will form a large part of the future measurement pro- 
gramme. Our experimental technique is ultimately limited by the 
precision of the control and measurement of the magnetic fields in 
the atom trap and its surroundings. Offline magnetometry using elec- 
trons, nuclear magnetic resonance** (NMR) probes, and possibly 
trapped, laser cooled ions‘, will lead to refinement of the current 
method. The central trapping region of ALPHA-g, not yet utilized, is 
designed to be less susceptible to unprogrammed magnetic fields 
and to work with colder atoms. Having a cold source of stable antimat- 
ter in a vertical trap suggests the possibility of performing fountain- 
type, gravitational interferometry measurements”, promising preci- 
sions of order 10% inthe determination of ag. Formerly the subject of 
countless thought experiments and indirect inferences, the motion 
of antimatter in the gravitational field of the Earth finally has a sound 
and promising experimental foothold 
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Methods 


Detection of antihydrogen annihilations 

The ALPHA-g radial time projection chamber. The rTPC is a three- 
dimensional (3D) particle tracking detector, designed to reconstruct 
the antihydrogen annihilation location from the charged m-mesons 
released in the process”. The detector has a cylindrical structure, 
placed between the ALPHA-g trap and the1 T solenoid magnet (Fig. 1). 
The ionization electrons created as a charged n-meson passes through 
the gas medium (mixture of 70% Ar and 30% CO.) drift to the outer 
walls of the detector, where they are collected, generating a read-out 
signal. The axial (z), azimuthal (@) and radial (r) position information 
about the particle trajectories is inferred from the signals induced on 
the segmented cathode pads (4 mm z pitch) and anode wires (4.5 mm, 
or 1.4° @ pitch), as well as from the drift time—the time it took the elec- 
trons to reach the outer wall (typically on order of microseconds). The 
1.8 x 10°cm? active gas volume is 230 cm tall inz, and extends from the 
inner cathode wall (r= 10.9 cm) to the segmented outer cathode wall 
(r=19.0 cm). The gas volume consists of two regions: a drift region 
(r=10.9 to 17.4 cm), where the main tracking information is obtained, 
and a proportional region (r= 17.4 to 19.0 cm), where electron multipli- 
cation takes place, inducing signals on 256 anode ‘sensing wires’ and 
onthe outer cathode pads. The pads have a 576-fold segmentation in z 
and 32-fold in @ (11.25°), for a total of 18,432 readout channels. A radial 
drift electric field (E,) is applied orthogonal to the axial solenoidal 
magnetic field (B,), making this a relatively uncommon configuration 
for a TPC**“’. This design choice was driven by factors including (1) the 
large aspect ratio of the height (approximately 230 cm) to the radial 
width (approximately 10 cm) of the available space, (2) the influence 
of the non-uniform magnetic fields (from the internal magnets and 
the solenoid fringe field) on the charge drift and (3) the capability to 
operate the detector at a lower or zero B,. Due to this field configura- 
tion, an electron that drifts radially outwards due to £, is also displaced 
in ġ, when B, is present (Lorentz displacement). The angular deflection 
is around 9° for maximal drift length at B, =1T. 

The barrel scintillator. The barrel scintillator (BSC) surrounds the 
rTPC and provides additional information on annihilation events. In this 
work, it was mainly used to provide information on the event topology, 
as a part of the cosmic background rejection analysis (see below). The 
BSCis composed of 64 trapezoidal scintillator bars (Eljen Technology 
EJ-200) that are 2.6 m long and 2 cmthick. The bars are read out at each 
end by anarray of six silicon photomultipliers of type SensL J-series, 
each photomultiplier having an active area of 6 mm x 6 mm. The ana- 
logue signals from the six SiPMs at each end (top or bottom) of a BSC 
bar are summed ina front-end card on the detector and sent through 
5 mcoaxial cables to a digitizer module and atime-to-digital converter 
for each of the 128 channels. 

Reconstruction. A charged t-meson typically produces about three 
ionization clusters per millimetre of track length in the rTPC drift 
volume. The determination of the 3D position (space point) of the 
cluster from the detected signals in the pads and the wires requires 
a model of the charge drift process in the detector gas medium. We 
use a sSimulation™ based on Garfield++ that accounts for Lorentz 
displacement in the non-uniform £, and B, fields for a given gas con- 
dition. In our rTPC configuration, multiple clusters from the same 
incoming particle often register signals on the same wires or pads. 
Atemplated-based deconvolution method was used to infer the space 
points for these events. Given the set of space points, the particle 
trajectories are identified using an algorithm that finds the nearest 
neighbouring point. A least-squares method is employed to find the 
best fit with a functional form of a helix in three dimensions. Finally, 
the antihydrogen annihilation position, or vertex, is calculated by 
finding the point where at least two helices pass closest to each other”. 
Extended Data Fig. 4 shows a distribution of the reconstructed ver- 
tices from antiproton annihilations in a Penning trap, indicating az 


vertex resolution of order 2 cm, which is considerably smaller than 
the 25.6 cm separation of the mirror coils A and G. 

Machine-learning suppression of the cosmic ray background. 
Cosmic rays are the dominant source of background. The cosmic event 
rate of around 70 Hz is suppressed by a factor of approximately 350 
by offline machine-learning analysis. Twenty selection variables that 
are sensitive to the topological differences between annihilation and 
background events are used as inputs to a boosted decision tree classi- 
fier". The machine-learning classifier is trained using experimental 
datasets of signal and background events. The signal sample (371,362 
events) was obtained from antihydrogen produced during antiproton 
and positron mixing near z = 0, filtered to only include periods of high 
event rates to minimize cosmic contamination. The background sam- 
ple (610,942 events) was collected when there were no antiprotonsin 
the apparatus. None of the machine-learning variables in the training 
samples showed significant correlation with the vertical axis (z). A 
classifier cut was chosen to optimize the significance for an expected 
1count of signal in a period of 8 seconds. This yields a detection effi- 
ciency of 0.65 + 0.02 annihilations per readout trigger. The background 
rates depend slightly on z and are reported in the Table 1 caption. 

Top and bottom detection efficiencies. The gravity measure- 
ment requires knowledge of the relative detection efficiencies for 
antihydrogen released in the defined ‘up’ and ‘down’ regions of the 
trap. The t-mesons from antihydrogen annihilating in these regions 
traverse slightly different amounts of material and magnetic fields. As 
discussed in the main text, the +10g datasets provide a direct calibra- 
tion of the relative detection efficiencies, because the vast majority of 
annihilations occurs in only one region. We have performed two other 
calibration measurements: (1) similar to the above but with a slightly 
different antihydrogen release sequence and (2) ameasurement where 
only one mirror was ramped downat atime, with the other being held 
at full field. These measurements gave a consistent relative efficiency, 
corroborating our +10g measurement. 


Detector performance and laser calibration. The stability of the 
detector high voltage, gas flow, temperature, and pressure were closely 
monitored during the measurement campaign; no trends that would 
affect the detector efficiency were observed. Regular accumulations 
of cosmic ray events were taken to monitor detector occupancy, noise 
levels and background rates. Throughout the campaign, 100% of the 
rTPC wires, more than 97% of the rTPC pads and more than 99% of 
the BSC channels were active. Data from faulty or noisy channels 
were removed from the analysis. 

A dedicated calibration system was developed to validate the Gar- 
field++ charge drift simulation. A266 nm pulsed laser illuminated nine 
aluminium strips (6 mm wide) placed along the inner cathode of the 
detector. This generates photoelectrons at well-defined z and @ posi- 
tions and at known times. Extended Data Fig. 5 shows good agreement 
between the calibration data and the simulation. The calibration also 
served to monitor variations of drift time influenced by environmental 
conditions throughout the measurement campaign. The track recon- 
struction analyses, performed by artificially varying the Garfield++ 
model values within the range indicated in Extended Data Fig. 5 (top 
and bottom panels), did not produce any significant changes in the 
reconstructed vertices, confirming the validity of our understanding 
of the detector and its robustness against the possible variations in 
operational conditions. 


Field measurement and modelling, magnetic biases 

Electron cyclotron resonance. In electron cyclotron resonance (ECR) 
magnetometry the magnetic field is deduced from the response of a 
test cloud of electrons to microwave radiation near its cyclotron fre- 
quency. The temperature of a single such test cloud, subjected toa 
single frequency of microwave radiation, is destructively measured 
through slow extraction to a microchannel plate and phosphor screen 
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assembly”. A spectrum can then be mapped out by rapid repetition of 
such single exposures using the reservoir technique while sweeping 
the microwave frequency. We fit a Gaussian function to the spectrum 
to extract the peak frequency. Here, we apply no evaporative cooling 
to the test clouds before exposing them to microwaves. This serves to 
minimize the radial extent (around 0.1 mm) of the test clouds and con- 
sequently their sensitivity to radial field gradients. This is necessary in 
the highly inhomogeneous magnetic fields in regions of the trap that are 
crucial to the current work. The microwave radiation is produced bya 
Keysight E8257D synthesizer, with a frequency resolution of 0.01 Hzand 
an amplitude accuracy, in the parameter range of interest, of +1.3 dB. 


ECR during the measurement trials. During each of our experimen- 
tal trials, we measure the magnetic field simultaneously at two fixed 
locations near the axial centres of the A and Gcoils immediately before 
and after the ramp-down of the mirrors, and again after zeroing all 
the currents in the internal magnets. This last measurement serves to 
monitor the stability of the background field over the course of many 
such trials. The simultaneous measurements are achieved by extract- 
ing and positioning two test clouds at atime and irradiating both with 
the same microwave pulse. 

The measurements before the mirrors ramp down display a broad- 
ened spectrum due to the high field gradient and have a full width at 
half maximum of order 7 x 10°T, while measurements in the final well 
and background (external solenoid only) magnet configurations havea 
full width at half maximum of order 2 x 10°T. While significantly smaller 
spectral widths can be achieved by tuning the microwave parameters 
and the test clouds, the settings used in this work were chosen to encom- 
pass many of the current configurations in the same linear frequency 
sweep and to ensure robustness against small changes in the loaded 
reservoir across many experimental trials. 


Rapid cycle ECR measurements. The repetition rate of obtaining an 
ECR spectrum is limited by the time it takes to load and prepare the 
reservoir from which we extract the test clouds. For the measurements 
before and after the gravitational release ramps, we extract, expose, 
and dump 200 clouds to produce simultaneous field measurements 
near the mirror A/G coil centres in 67 s. By comparison, the reservoir 
was loaded in around 75 s. A faster repetition rate was obtained by using 
areduced set of 25 microwave exposure frequencies to produce eight 
repeated measurements from the 200 test clouds. In addition, with 
careful tuning of the reservoir and the test cloud extraction, we can 
also extract more test clouds froma single reservoir; see the magnetron 
frequency magnetometry section below. As illustrated in Extended Data 
Fig. 6, we used this technique to track the decaying field immediately 
after the end of ramping down the mirror coils. The reservoir was loaded 
during the magnet ramp, so the resonance was hit within 3 s of the ramp 
completing. The fits to the data are sums of two exponential decays 
with differing time constants (roughly 20 s and 300s). 


Field measurement using the electron magnetron frequency. We have 
developed a technique that uses the magnetron frequency ofan electron 
plasma as a measure of the magnetic field at various axial positions inthe 
ALPHA-g device. The measurements described below are taken offline. 
Using the reservoir technique‘ we extract two thousand reproducible 
‘electronclouds’, each containing about 1,000 electrons at a temperature 
of 100 K and a radius of 100 um. Although patch potentials (unpro- 
grammed potentials due to, for example, charged oxide layers) and volt- 
age offsets cause the trapping potential to differ from an electrostatic 
model by about 1%, these potentials are reproducible from day to day to 
at least one part in 10°. When a cloud is radially displaced from the trap 
centre and trapped by an electrostatic potential V+ approximated 
byV; (z, r) =k,(z? - r?/2), where k, is determined by the electrode poten- 
tials, it orbits the centre at a frequency w, , given by w, = k,/|B|. Precise 
measurements of this frequency are performed in the following way: 


1. Acloud is extracted from the reservoir and moved axially to the 
desired measurement location. 

2. Patch potentials introduce a transverse electric field in the otherwise 
cylindrically symmetric Penning trap. When the trapping potential 
is weak (O[0.5 V]), the magnetron orbit is no longer centred****. We 
quickly decrease the trapping voltage and wait (about 10 ms) for the 
cloud to arrive at a desired off-axis location. 

3. The trapping voltage is then quickly increased, and the cloud begins 
to orbit the trap centre. After a variable amount of time, it is released 
towards one of the multichannel plate (MCP) detectors (Fig. 1). The 
final magnetron phase can be extracted from the cloud’s imaged 
position. 

Asingle image does not suffice as ameasurement of the magnetron 
frequency because the cloud’s total number of orbits is ambiguous. 
First, we image one cloud that orbits the trap fora short time T,= 100 us. 
Then we image N clouds that orbit in the 1 T magnetic field for atime 
T,= Ty + 1.4°1(1T)/k,, for n=1 to N, that is, geometrically increasing 
the hold time. For several reasons, there is a variability in the final 
angular position of about 0.1-0.4 radians depending on the axial loca- 
tion of the measurement. The constant 1.4 is chosen such that before 
each measurement, the estimate of w, is good enough that there will 
be no ambiguity in how many times the cloud orbited the trap centre. 
In this way, we can increase the total magnetron phase angle while 
having a roughly constant error. To extract the magnetic field from 
a precise measurement of w,, we calibrate w, at a particular field 
measured with ECR, and we use the relationship œ, « 1/|B| to measure 
the field inthe presence of different magnet currents. Of course, there 
are corrections to this relationship, which are at most about one part 
in 10*. 

This technique has been useful for measuring the magnetic field 
while the magnet currents are ramping. To do this, we image succes- 
sive clouds after a time T,= 2,000(1 T)/k, (that is, an amount of time 
such that the cloud would orbit about 2,000 radians in a1 T field). As 
the field decreases, the magnetron frequency increases. Depending 
on the location in the trap, we perform measurements once every 
30-50 ms, which means that, in a 20 s ramp, each cloud orbits the 
trap at most 5 radians more than the previous cloud. We track the 
total magnetron angle by initially employing a ‘geometrical increase’ 
operation before the field changes, then we add the angle deviation 
between successive clouds. 

Extended Data Fig. 7 shows an example measurement of a 20 s mag- 
net ramp inthe centre of mirror A. The first subplot shows the raw 
measurements. The second shows the result of subtracting the 
expected model for the magnetic field, which assumes it changes lin- 
early between ECR-measured magnitudes before and after the ramp. 
The most striking feature is a nonlinear component of about 1 x 10°T, 
which we interpret as persistent currents being induced into super- 
conducting material. When a magnet’s current is decreased by A/ from 
a starting value /,, we observe a nonlinear component of the field that 
exponentially saturates with increasing A/. For the mirror coils 
ramp-down used to measure gravity, the field is well approximated by 
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By performing this magnetron magnetic field measurement during 
130 s magnet ramps in 20 axial locations, we measured A(z), and this 
behaviour of the magnetic field was added to antihydrogen simula- 
tions (see below). The A(z) produced by a mirror coil ramp-down looks 
similar to the nominal field produced by the mirror coils; it has two 
bumps centred on each coil. In other words, persistent currents resist 
a change in the magnetic field. We measure a small difference in A(z) 
at the locations of mirror coils A and G that gives rise to the approxi- 
mately exponentially saturating behaviour of the bias at early times in 
Extended Data Fig. 6 (see below). 


The final subplot of Extended Data Fig. 7 includes three corrections: 
1. The exponentially saturating persistent field used in the magnetic 

field model is subtracted. 

2. Toimage off-centre clouds onthe MCP, additional normal conducting 
magnets near the MCP need to be energized. The 0.6-1.0 mT effect 
of these magnets is subtracted. 

3. The frequency w, depends on the distance a cloud is displaced from 
the trap centre—in part because |B| increases off-axis. The correction 
(about five parts in 10°) from this effect is obtained from separate 
calibration measurements. It takes the form mA/ + b, so the final 
subplot includes constant, exponentially saturating and linear 
corrections. 

Despite these corrections, the field shows some deviation from the 
‘expectation’. First, the deviation is about -0.1 mT before and after 
the field starts changing. The most likely explanations are errors 
in the measurement technique that are linear in|B| (including calibra- 
tion error). While the field is changing, there is a positive deviation of 
0.1 mT. This is a known effect from the induced current in a nearby 
magnet. Next, there are exceptional measurement points just after the 
magnet ramp starts and just after it ends; these are known effects of 
the magnet control system. There is also a small increase in the first 
second because the persistent current is not perfectly modelled by an 
exponentially saturating function. Only this last effect is not included 
in simulations of the experiment, but it occurs inthe same way in both 
mirror coils and so does not affect the bias. Inthe end, the magnetron 
technique provides certainty that there are no other unmodelled 
effects in the on-axis magnetic field larger than 0.1 mT. 

Similar data were taken for several biases at five locations near the 
centre of each mirror coil. Additionally, the magnetron technique 
was used to measure magnetic fields in 20 axial locations through- 
out the trap during the 130 s magnet ramp-downs. These data were 
useful for identifying and quantifying the exponential saturation of 
persistent currents. The longer measurement time allowed for amore 
precise measurement of A(z), which we later verified was consistent 
with what we observe in 20 s ramp-downs. An upcoming publication 
will provide a more detailed analysis of these data and description of 
the measurement. 


Bias uncertainties. Table 2 lists the estimated uncertainties in our 
calculation of the on-axis bias. Here we detail how each of those contri- 
butions is estimated. Firstly, each ECR spectrum taken exhibits a finite 
width constituting an uncertainty in the determination of the magnetic 
field from that spectrum”. Since the magnetic field difference (B, - By) 
is what enters the bias, we add in quadrature the fitted Gaussian widths 
from measurements in mirrors A and G. We then average over all valid 
ECR measurements at the beginning and end of the release ramp to get 
the ‘ECR spectrum width’ contribution. 

The ‘repeatability of (B, — Bay contribution describes how well the 
magnetic field difference is repeated from one experimental trial to the 
next and is evaluated as the standard deviation of all valid bias measure- 
ments around the average in each set. 

Due to background field gradients caused mainly by the octupole 
windings, the on-axis field maxima at the end of the ramp are shifted 
away from the geometric centres of the mirror coils as the currents 
decrease. We correct for this by mapping out the field maxima with 
high spatial resolution for every current configuration used (Fig. 2a). 
Parabolic fits are then used to extract the true locations of the on-axis 
field maxima (saddle points in 3D), as well as the difference between 
the field measured at the two fixed locations during the gravity experi- 
ments and the true maxima. We take the average absolute residuals of 
the parabolic fits as an error in this correction, adding in quadrature 
the errors evaluated in the two mirrors and averaging over all current 
configurations. This is tabulated as ‘peak field size and z-location fit’. 

The ‘field decay asymmetry (A to G) after ramp’ uncertainty arises 
because there is a delay (about 96 s) between the end of the mirror ramp 


and the measurement of the magnetic field. We expect a slight change 
in magnetic field in this time due to the decay of persistent currents 
induced by ramping the magnets. If this decay is not equal in the two 
mirror coils, there would be an error in the field difference measured. 
The fast repeat ECR described above allowed us to quantify the field 
decay and look for any asymmetry in a dedicated measurement that 
is shown in Extended Data Fig. 6. Here we shift the data to overlap the 
fitted fields at O s and to best highlight any difference in decay rate. 
We observe a 6 x 10°°T field change during the first 96 s after stopping 
the ramp, with no appreciable asymmetry between the two mirrors, 
nor a strong dependence on the exact current configuration. We take 
as a potential error the largest observed decay difference between 
the mirrors out of the three biases investigated. 

In the main text, we describe how the time-averaged bias for each 
current configuration is calculated by averaging the calculated bias 
present in the trap at the time of each annihilation event. This is illus- 
trated in more detail in Extended Data Fig. 8, for anominal bias of Og. 
The uncertainty we associate with this spread of biases is the standard 
deviation of the individual calculated biases of annihilation events. The 
number given for the ‘bias variation in time’ uncertainty in Table 2 is 
averaged over all current configurations; individually, they range from 
0.010g to 0.035g. 

The bias calculations above rely ona field model to extract the bias at 
any time during the ramp. The field model is constrained both by ECR 
measurements of the field at various currents as well as magnetron 
frequency measurements (see below). To evaluate the accuracy of the 
on-axis bias inthe model, we compared it to offline (that is, independ- 
ent of the experimental gravity trials) ECR measurements taken in both 
mirrors at 10 points along the nominal magnet ramps, making sure 
to match the magnet ramp history and resulting induced persistent 
currents to the gravity trials. We repeated these measurements for 
five different current configurations and define the global average of 
absolute residuals to be the ‘field modelling’ uncertainty. 


Simulations of the dynamics of trapped antihydrogen 

Field model. A field model was developed to include all know- 
ledge of the magnetic trap during the mirror A/G ramp-down. The 
model was used to derive the on-axis trap biases and to simulate the 
three-dimensional trajectories of atoms in the trap. 

For the external (1 T) solenoid, an ideal field was first calculated from 
the designed winding geometry. This was compared to field measure- 
ments made with a rubber sample NMR probe in the empty solenoid 
bore. The difference between the two was deconvolved, using singular 
value decomposition, to yield current density perturbations on the 
solenoidal windings. The subsequent installation of the inner cryostat 
and coils into the external solenoid perturbed its field. The change, 
mapped on-axis by ECR, was deconvolved into a model solenoidal cur- 
rent distribution overlapping the inner superconducting windings. The 
ECR-measured background field was replicated in the field model to 
within 5 x 10°T. Inthe simulation, this background field was assumed 
to be static during the A/G ramp-down. 

The field contributions from the octupoles and mirror coils were 
computed from winding geometries measured during fabrication. 
The model windings were slightly offset and scaled to best match the 
ECR mapping of individual magnets. The currents used in the field 
model during the A/G ramp-down were measured experimentally 
using direct-current current-transformers (DCCTs). The experimen- 
tal current histories had a sample rate of 10 kHz and were filtered by 
removing Fourier components above 1 kHz before being applied in 
the field model. 

Field measurements made during the mirror A/G ramp, with all 
windings energized together, revealed field contributions that did 
not originate from the applied current in any individual winding. We 
model these contributions in two parts: an exponentially saturating 
component derived from the magnetron measurements described 
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above, and a residual linear component that further improves the agree- 
ment with the aggregate field measurements. These contributions, 
approximately 10°T in magnitude, were included in the field model 
usinga time- and z-dependent solenoidal current distribution located 
approximately where the inner superconducting windings are located. 

Putting all contributions together, the field model produced fields 
that agreed with online ECR, offline ECR and magnetron measurements 
toa standard deviation of around 2 x 10°T overall and around 1 x 10°°T 
near the trap saddle points at coils A and G. The former value, converted 
to units of bias, is quoted in Table 2 as the ‘field modelling’ uncertainty. 


Trajectory computation. The field magnitudes were precomputed, 
stored in a regular grid of 0.5 mm spacing, and interpolated via a 
third-order polynomial for the trajectory simulation. The field inter- 
polation was fractionally accurate to 10% near the cylindrical vacuum 
wall where the fields had high spatial variations, and was substantially 
better away from the wall. 

Atoms were evolved in time using a leapfrog stepping algorithm. The 
time step was chosen individually for each atom and was either 1 us 
or an interval such that length traversed per step was no longer than 
0.03 mm at all times, whichever was smaller. Stepping was terminated 
when atrajectory reached the inner Penning trap electrode surface, the 
UHV beam pipe, or two artificial axial stops located outside the region 
where atom annihilations are registered by the detector. 

The trajectory simulation was made in two parts. (1) To model the 
initial catching and accumulation process, atoms were initialized near 
the bottom of the trap. The positions were uniformly distributed over a 
cylinder of 1 mm radius and 5 mm length. The velocity was drawn froma 
50 K Maxwellian distribution. The atoms were initialized with a principal 
quantum number of 30 and allowed to radiatively cascade down to the 
ground state using the method described by Topçu and Robicheaux”. 
Each atom was evolved for a randomly selected duration between 0 and 
14,400 sto simulate the gradual accumulation of antihydrogen during 
‘stacking’. The 6,726 atoms that remained trapped after their specified 
duration were retained. (2) These atoms were evolved in time through 
the long octupole and the A/G mirror coils ramp-down using various 
trap biases and under various assumed gravitational accelerations. 
The time and location of annihilation were recorded, from which the 
escape bias curves in Fig. 5, Extended Data Fig. 1 and Extended Data 
Fig. 2 were derived. 


Systematic uncertainties. In addition to the escape curves, other 

results from the simulation have been compared to the experiment. 

The escape time and axial position distributions of annihilation verti- 

ces during the LOc and mirror A/G ramp-down windows showed good 

agreement. On the other hand, the behaviour of atoms that remained 
after the A/G ramp-down differed. (Note that these atoms do not con- 
tribute to escape curves.) Inthe simulation, one annihilation in the LOc 
window corresponded to 0.08 annihilations during the hold after the 

A/G ramp-down, and 0.51 during the subsequent OcB ramp-down. In 

the experiment, these numbers were 0.27 and 0.10. This meant fewer 

atoms than expected survived the A/G ramp-down, and more atoms 
were driven out of the trap during the hold despite the trap field remain- 
ing nominally unchanged. 

Given these differences, parameters in the simulation were per- 
turbed to establish the robustness of the escape curve, and to obtain 
the uncertainty shown in Fig. 5 and quoted in the measured value of 
the antihydrogen acceleration towards the Earth. We considered the 
following: 

1. The disagreement in the fraction of atoms surviving the A/G ramp- 
down was found to be consistent with the simulation not having 
initialized the atoms’ energy in the same way as the experiment. As 
in our previous work”, uniform and linear initial energy distributions 
were simulated by bootstrapping the results of the nominal 50 K 
Maxwellian initial energy simulation. The escape curves resulting 


from these distributions tended to have lower central slopes com- 
pared to the nominal curve, but the point of balanced escape re- 
mained unchanged. The uncertainty in the simulated escape curve 
due to this analysis of the total initial energy distribution is included 
inthe uncertainty band in Fig. 5. This demonstrated that the escape 
curve was not sensitive to even drastic changes to the initial condi- 
tion of the atoms. 

2. The higher-than-expected annihilation count during the hold after 
the A/G ramp-down was consistent with an energy exchange between 
the transverse and parallel degrees of freedom that was not pre- 
dicted. Anartificial, unphysical exchange mechanism was therefore 
introduced to the simulation where atoms received random velocity 
deflections during their evolution. The strength of this artificial 
deflection was constrained by the timing of escapes, as excessive 
exchange forced atoms to escape early. Within this constraint, no 
changes to the escape curve were observed. 

3. Multipolar perturbations with zero component on axis can alter the 
escape balance of the experiment while eluding ECR and magnetron 
measurements. Dipole, quadrupole, sextupole and octupole field 
perturbations were applied to the bottom half (z < 0) of the trap to 
maximize the induced asymmetry. Assuming these perturbations 
arose fromerror inthe radial positioning of the OcB conductor, the 
multipolar fields were constrained by the accuracy with which the 
winding was fabricated (around 10 um). Assuming the field pertur- 
bation arose from persistence effect, the multipolar fields were 
constrained by the critical current of NbTi. The former resulted ina 
stronger perturbation and was simulated. The octupole mode per- 
turbation had the most significant impact on the escape curve and 
effected a maximum +0.26g offset along the bias axis. The central 
slope was unchanged by the perturbations. The uncertainty (one 
standard deviation of an assumed flat distribution, Table 3) in the 
simulated escape curve due to the octupole mode perturbation is 
included in the orange uncertainty band in Fig. 5. 

4. Other field perturbations that were consistent with on-axis mag- 
netometry measurements included transverse offset of the axis of 
the A and G coils from the OcB axis, and angular misalignment of the 
external solenoid. These resulted in no change to the escape curve 
within the mechanical constraints. 

5. Mechanical vibration of the trap magnets could heat the trapped 
atoms and alter their dynamics. This was simulated and no changes 
to the escape curve were observed at vibration amplitudes below 
obviously audible/tactile limits. 

For each bias value on the escape curve, the largest positive and 
negative deviations from the unperturbed P,, resulting from the above 

perturbations were chosen for the band displayed in Fig. 5. 


Magnets and magnet controls. The ALPHA-g magnetic trap is gener- 
ated from superconducting windings housed in two cryostats: the outer 
cryostat houses a solenoid and shim coils that provide the uniform axial 
background magnetic field of 1 T needed for plasma confinement inthe 
Penning trap, while the inner one contains 21 distinct superconducting 
circuits ©°. Figure 1 in the main text shows the subset of magnets in use 
for this study. Mirrors A and G are used to provide axial confinement 
to the anti-atoms as well as to control the release and are energized in 
series up to approximately 70 A by a CAENELS FAST-PS-1K5 operating 
in voltage controlled current supply mode (16-bit analogue to digital 
input with analogue bandwidth of 1 kHz). An additional, much smaller, 
differential currentis applied in parallel to mirror G alone, using a Kepco 
BOP 20-10 in voltage controlled current supply mode (analogue input 
with 10 kHz bandwidth) (Extended Data Fig. 9). We label the series and 
differential circuits as MAG and MGDiff respectively. This connection 
scheme ensures that any noise or drift in MAG is shared between both 
coils and thus has a small impact on the up-down balance of the trap. 
Extended Data Table 1 details the power supply and performance char- 
acteristics of the circuits used in the atom trap region. 


We use PM Special Measuring Systems TOPACC Zero-Flux DCCTs 
installed on the magnet current leads to actively monitor the current 
supplied to the magnets. The MGDiff circuit was measured using 30 
turns of its lead through its DCCT head. Calibrated accuracy of the 
units is about 25 ppm of the DCCT’s full scale (around 2.5 mA-turn), 
with less than 1 ppm drift expected over the course of this experiment. 
Full-scale output of the DCCT is transmitted by a +10 V signal with an 
output small-signal bandwidth of 500 kHz. The DCCT output voltages 
were digitized with 24-bit +10 V National Instruments NI-9239 cRIO ADC 
modules ata rate of 50 kS s”. Firmware on the NI cRIO FPGA recorded 
arunning average of this signal ata rate of 10 kS s”. This measurement 
was used for proportional-integral—derivative (PID)-based closed-loop 
control of the magnet power supplies (excluding the external solenoid 
supplies) to compensate for non-linearities in the QPU circuits and 
internal drift of the power supplies. Current programming voltages for 
power supplies were generated by NI-9264 analogue output modules 
with 16-bit resolution. Parallel readout of all monitored and control 
voltages was recorded at 10 kS s by the firmware, with jitter on the 
order of 1 ps and clock drift relative to the main data acquisition system 
at the 10 ppm level. 

Currents measured during 20 s and 130 s ramp-downs achieved 
run-to-run repeatability within the operating noise level of the mag- 
net systems (Extended Data Table 1). Deviations from the requested 
current included a consistent and constant current offset of 1.5 mA 
during the 20 s linear ramps and 0.22 mA for 130 s linear ramps. These 
offsets were due to lag in the PID control loop. In addition, a consistent 
overshoot transient at the start and end of the ramps was produced by 
the PID control of the MAG series circuit. The deviations of the MAG 
series current from the programmed linear ramp directly affect the 
atom trap depth and also introduce a bias field error due to the roughly 
1% construction difference between coils. For the 20 s ramp, this was a 
swing of less than 80 mA (bias 0.017g) over approximately 200 ms at 
the start of the ramp and less than 65 mA at the end of the ramp (bias 
0.014g, or 12% of final well depth). For the 130 s ramp the start transient 
was less than 15 mA (bias 0.0032g) over 200 ms and less than 12 mA (bias 
0.0025g, 2% of final well depth) over approximately 200 ms. 

During release measurements, currents were inductively coupled 
into mirrors B, C, D, E and F (unpowered and shorted through resistors), 
though not inthe Background and Transfer coils (disconnected during 
this study). The respective currents in mirrors B through F were meas- 
ured during release ramp-downs to be 7.9 mA, 2.6 mA, 2.1 mA, 2.6 mA 
and 8.1 mA during 20 s ramps, and 1.2 mA, 0.4 mA, 0.3 mA, 0.4 mA and 
1.3 mA during 130 s ramps. These contribute to bias magnetic field 
errors at a level well below 0.01g. All measured currents were included 
in the numerical simulations of the experiment. 


Analysis for escape curve and gravitational acceleration. The analy- 
sis begins by aggregating the time and axial location of antihydrogen 
annihilations reconstructed during the mirror ramp-down for each 
bias. Next we apply the z and time cuts, described in the main text, to 
the data. Using experimental calibration samples with biases of -10g 
and +10g, for which antihydrogen is largely forced to escape upwards 
or downwards, we calibrate the efficiencies in the up and down regions 
of the detector. The cosmic background rates across the trap are con- 
strained using data obtained while the trap is empty. 

We perform a likelihood analysis“ to determine the probability to 
escape downwards, P,,, (or equivalently the asymmetry A between the 
downward and upward escaping anti-atoms A = 2P,,,— 1), at each bias. 
The credible intervals for P,,, are shown in Fig. 5. 

Using the simulation, we then find the set of simulated downward 
escape probabilities, P,;,,, at the measurement biases, for a range of 
simulated values of the gravitational acceleration asim and perform a 
likelihood analysis on the experimental data to estimatea z. The results 
are cross-checked by repeating the analysis with different fiducial cuts 
in tand zand with the 130 s ramp data. 


We estimate the significance of having observed the effect of gravity 
on antihydrogen from the asymptotic distribution of the likelihood 


ratio between the models with zero and the extracted value of az. 


Counting statistics are included in the likelihood analysis by assum- 
ing that the counts in the mirror release in the up and down regions and 
the LOc counts at each bias are sampled from independent Poisson 
distributions with the mean specified in terms of the experimental 
parameters. 

Systematic uncertainties are included by allowing the parameters 
that enter the likelihood analysis to vary according to their experimen- 
tal uncertainties (where available) or within plausible ranges. The 
dominant source of systematic uncertainty in estimating P,, is the 
calibration of the detector efficiencies in the up and down regions. 
The dominant source of error in calculating az is related to errors 
inthe simulation model arising from uncertainties in the off-axis mag- 
netic field. Table 3 provides a breakdown of the contributions consid- 
ered for the total uncertainty. 


Data availability 
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Extended Data Fig. 1 | Simulated escape curves for various values ofa z. We for details. The simulations are otherwise identical to that used for the 20 s 


illustrate the escape curves resulting from assuming several different valuesof release experiment and normal gravity (Fig. 5). The solid blue curve represents 
the gravitational acceleration, az, of antimatter due tothe Earth. Seethelegend expectations for normal gravity. 
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Extended Data Fig. 2 | Escape curve for 130 sramp-down. The P,, values are those for the 20 s protocol. The 20 s data and simulations for both ramp times 
plotted versus bias for the three trial sets having biases 0 g,-1g,and-2g.These are also shown for comparison. Note that the simulated escape curve for the 
biases were chosen after the 20 s ramp results had been examined. Apart from 130s ramp has a steeper transition region than for the 20 s ramp, and the 


the slower ramp, the experimental and analysis procedures were identical to balance point (Pan = 0.5) is not at a bias of precisely -1 g, as described in the text. 
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points are assumed to be the same as for the blue points. This assumption 
should be valid within the uncertainties. The vertical error bars are obtained 


by following the same procedure described in the text. The blue (green) curve 
is based on 1,722 (621) total events as defined in Table 1. 


Extended Data Fig. 3 | Escape curve for atoms escaping after 20 s.Theescape 
curve (green points) for the time period (20-80 s) after the mirrors AandG 
have stopped ramping down and are held at constant current while the ECR 
measurement is prepared. The main (10-20 s) data set (blue points) is shown 
for comparison. Note that the bias values and their uncertainties for the green 
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Extended Data Fig. 4 | rT PC resolution: a measured z-distribution of with two Gaussians and a flat background. The two distributions have 
annihilation vertices from antiprotons held ina Penning trap. An standard deviation 1.5 cm (Gaussian 1; -70% of the counts) and 4.2 cm (Gaussian 
approximately point-like source of events is obtained from antiprotons 2,~-24% of the counts). Both widths are significantly smaller than the distance 
annihilating on residual gas while being held for 2000 sinashort Penning trap. between mirrors A and G (25.6 cm, magnet centres illustrated by the green 


The reconstructed vertex distribution in z (points with error bars) is fitted vertical lines). 
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Extended Data Fig. 5| Laser based rT PC calibration and comparison wit 
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h dashed lines indicate the axial midpoints of mirrors A and G. b The red circles 


Garfield++ simulation predictions. a The green bars denote the measured denote measurements, at five z-locations, of the Lorentz displacement, i.e., the 


drift times of laser-induced photoelectrons released from nine aluminium 


azimuthal displacement of the radially drifted photoelectrons at the location 


strips onthe inner cathode surface of the rl PC (Methods). The green dashed of the anode wires (r=18.2 cm). The error bars are due to the fit error and are 


line is the predicted drift time from the Garfield++ model used inthe detec 
physics analysis. Each green bar denotes the range of drift times measured 
during the period of the physics measurement campaign, while the dashin 
each bar denotes the average drift time over those measurements. Vertical 


tor smaller than the plotting symbols except at the vertical position O cm. The two 
points for each zare measured at different initial azimuthal positions. The red 
dashed line denotes the prediction from the same Garfield++ simulation 
model. Note that vertical scales are magnified for both figures. 
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Extended Data Fig. 6 | Decay of persistent fields (offline measurements). 
The on-axis fields at the axial midpoint (Fig. 2a) of mirrors A and G, as measured 
by the rapid cycle ECR technique (Methods) to study the decay of persistent 
fields after the end of the 20 s ramp. The solid lines are fits using two exponential 
decay times per curve; see Methods. The three plots are for the extreme biases 
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+3 gand Og. The red points represent the extracted systematic error inthe 
magnetic field difference between mirror A and mirror G at each bias, and they 
are plotted at the approximate time of the ECR resonance in the actual gravity 
trials. ‘Offline’ refers to measurements taken independently of the release 
experiments. 
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Extended Data Fig. 7 | Magnetic field measurements via the magnetron expected linear ramp is plotted versus time. c Measurements after accounting 
frequency. Magnetic field measured using the magnetron frequency of for the three corrections described in Methods. The 1-2 s gapsin the data are 
electron clouds inthe centre of mirror A versus time during the magnet dueto a memory limitation in the FPGA that controls the electrode voltages. 
ramp-down. a Raw measurements (blue) are compared to the expected linear New voltage instructions are loaded in this time. 


ramp (green line). b The difference between the measurements and the 
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Extended Data Fig. 8| Time dependence of the bias. Thesolidcurverepresents ofevents detected asa function of time for the O gtrials, asin Fig. 4. The red 


the modelled deviation of the bias from the nominal value (in this case Og) asa point shows the derived bias and its uncertainty. Note that the bias deviation 
function of time during the mirror ramp-down. Thehistogramshowsthenumber duringthe time data are collected is less than 0.1 g. See Methods. 
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Extended Data Fig. 9 | Schematic diagram of the circuit for energizing mirrors A and G (MAG) inseries and supplying the differential current (MGDiff) to 
mirror G only. The power supplies are described in Methods and Extended Data Table 1. QPU stands for ‘quench protection unit’. 


Article 


Extended Data Table 1| ALPHA-g magnet power supplies 


Supply Current Programming Current Bias field 
(maximum | Programming Resolution noise (mA | from current 
Operating ae (g) rms) noise 
Current 


Long Octupole Sorensen 37 TA 0.5 
(LOc) SGA 10- 
1200 


Bottom Octupole | Sorensen 
(OcB) SGA 10- 


CAENELS 
FAST-PS- 
1K5 

(70 A) 


Mirror G bias Kepco 
(MGDiff) BOP 20-10 
(3 A) 


External 2x 
Solenoid Main CAENELS 
Coil FAST-PS- 

1K5 

(191 A) 
External CAENELS 
Solenoid Shim FAST-PS 
Coil 1020-200 


The power supplies used for energizing ALPHA-g magnets are tabulated, together with some of their relevant performance characteristics. Resolution and noise values are expressed in units of 
the local acceleration of gravity for matter (9.81ms”). 


